BPS and non-BPS kinks in a massive non-linear S 2 -sigma model 
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The stability of the topological kinks of the non-linear S 2 -sigma model discovered in [l|] is discussed 
by means of a direct estimation of the spectra of the second-order fluctuation operators around 
topological kinks. The one-loop mass shifts caused by quantum fluctuations around the these kinks 
are computed using the Cahill-Comtet-Glauber formula 0. The (lack of) stability of the non- 
topological kinks is unveiled by application of the Morse index theorem. These kinks are identified 
as non-BPS states. There are two types of topological kinks coming from the twofold embedding of 
the sine-Gordon model in the massive non-linear sigma model. It is shown that sine-Gordon kinks 
of only one type satisfy first-order equations and are accordingly BPS classical solutions. Finally, 
the interplay between instability and supersymmetry is explored. 

PACS numbers: 11.10. Lm , 11.27. +d, 75.10. Pq 



I. INTRODUCTION 

The main theme in this paper is the analysis of the 
structure of the manifold of kink solitary waves discov- 
ered in In particular, we shall offer a full description 
of the stability of the different type of kinks. As a bonus, 
we shall gain information about the semi-classical behav- 
ior of such kinks from the stability analysis, providing us 
with enough data to compute the one-loop mass shifts 
for the topological kinks. 

Prior to our work [l[ , kinks in massive non- linear sigma 
models have been known for some time and profusely 
studied in different supersymmetric models under the 
circumstance that all the masses of the pseudo Nambu- 
Goldstone particles are equal. The study started with 
two papers by Abraham and Townsend [H in which 
the authors discovered a family of Q-kinks in a (1+1)- 
dimensional Af — (4, 4) supersymmetric non-linear sigma 
model with a hyper-Kahler Gibbons-Hawking instanton 
as the target space and mass terms obtained from di- 
mensional reduction. In however, these kinks were 
re-considered by constructing the dimensionally reduced 
supersymmetric model by means of the mathematically 
elegant technique of hyper-Kahler quotients. By doing 
this, the authors deal with massive CP N or MP N mod- 
els, a playground closer to our simpler massive § 2 -sigma 
model. Similar Af = 2 BPS walls in the CP 1 -model with 
twisted mass were described in [(| . In a parallel develop- 
ment in the (2+l)-dimcnsional version of these models, 
two-dimensional Q-lumps were discovered in Q and Q. 
Throughout this field, the most interesting result is the 
demonstration in Q and [Io| that composite solitons in 
d = 3 + 1 of Q-strings and domain walls are exact BPS 
solutions that preserve \ of the supersymmetries: ( Sec 



also the review llj, where a summary of these super- 



symmetric topological solitons is offered. 

Our investigation differs from previous work in the area 
of topological defects in non-linear sigma models in two 
important aspects: 1) We remain in a purely bosonic 
framework; in fact, we consider the simplest massive 
non-linear sigma model. 2) We study the case when the 



masses of the pseudo Nambu-Goldstone bosons are dif- 
ferent. The search for kinks in the (l+l)-dimensional 
model (domain walls in d = 3 + 1) is tantamount to the 
search for finite action trajectories in the repulsive Neu- 
mann system, a particle moving in an § 2 -spherc under 
the action of non-isotropic repulsive elastic forces. It is 
well known that this dynamical system is completely in- 
tegrable [25[, [I?! ■ We show, however, that the problem is 
Hamilton- Jacobi separable by using elliptic coordinates 
in the sphere. Use of this allows us to find four fami- 
lies of homoclinic trajectories starting and ending at one 
of the poles which are unstable points of the mechanical 
system. In the field-theoretical model the poles become 
ground states, whereas the homoclinic trajectories cor- 
respond to four families of non-topological kinks. Each 
member in a family is formed by a non-linear combina- 
tion of two basic topological kinks (of different type) with 
their centers located at any relative distance with respect 
each other. 

It is remarkable that the static field equations of this 
massive non-linear sigma model are (almost) the static 
Landau-Lifshitz equations governing the high spin and 
long wavelength limit of ID ferromagnetic materials. 
From this perspective, topological kinks can be inter- 
preted respectively as Bloch and Ising walls that form in- 
terfaces between ferromagnetic domains, similar to those 
discovered in the XY model dealt with in [33| . The vari- 
ety of our non-topological kinks, understood as solitary 
spin waves, is thus formed by non-linear superpositions of 
one basic Bloch wall and one basic Ising wall at different 
distances. Far from this non-relativistic context, degen- 
erate Bloch/Ising brancs have been studied in two-scalar 
field theories coupled to gravity in [l2|, [H, G3 ■ 



II. THE (l+l)-DIMENSIONAL MASSIVE 
NON-LINEAR S 2 -SIGMA MODEL 



We shall focus on the non-linear S 2 -sigma model stud- 
ied in Reference [l[ . The action governing the dynamics 
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with V = V(c/)i(t,x),(j>2(t,x),ip3(t,x)). The scalar fields 
are constrained to satisfy: 2 + <fi 2 + 0§ = R 2 , and thus 
4> a (t,x) G Maps(M ia ,§ 2 ) are maps from the (1 + 1)- 
dimensional Minkowski space-time to a S 2 -sphere of ra- 
dius R, which is the target manifold of the model. 

Our conventions for R 1,1 are as follows: x^ G R 1 ' 1 , 
H = 0, 1, x» ■ x„ = g^XpXv, g^ = diag(l, -1). x° = t, 



= t 2 -x 2 ; = g^d 2 ^ = □ = d 2 - d 2 
The infrared asymptotics of (1 + l)-dimensional scalar 
field theories forbids massless particles, see . We thus 
choose the simplest potential energy density that would 
be generated by quantum fluctuations giving mass to the 
fundamental quanta: 

V{4>1,4>2, h) = \ (a? 4>\ +al4>% + a\ 2 ) , (2) 

which we set with no loss of generality such that: a\ > 
a% > a§ > 0. 

1. Solving 03 in favor of 0i and 02, 03 = 
sgOfeV-R 2 -4>\ - 4%, we find: 



- 21^(01) 02 ) 



S = - / dtdx {d^id^i + d^ 2 d^ 2 - 

(01<9 M 0i + 02<9 / x02)(01<9 A '01 + 02<9 M 2 ) 



R 2 



V&a(0i,02) - 2 (("i ~ a i) 01 + ("a - "a) <^2 + const.) 



A 2 
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^iM + TT^M 



(3) 



with A 2 = (a 2 - a 2 ), 7 2 = («1 - a§), A 2 > 7 2 . 

2. Thus, the interactions come from the geometry: 

(01^0! + 2 d„0 2 ) {4>id^i + 02^0 2 ) 



i? 2 



■ (01^0! + 2 ^0 2 ) (01^0! + 2 C^0 2 ) 



and is a non-dimensional coupling constant, whereas 
the masses of the pseudo-Nambu-Goldstone bosons are 
respectively A and 7. 

Taking into account that in the natural system of units 
h = c = 1 the dimensions of fields, masses and coupling 
constants are [0 a ] = 1 = [R], [7] = M = [A], we define 
the non-dimensional space-time coordinates and masses 



a A 2 



to write the energy in terms of them: 

dx {(9 t 0!) 2 + (d t 2 ) 2 + (a,0!) 2 + (a,0 2 y 



(0ld t 01 + 02<9 t 2 ) 2 + (01(9x01 + 02 flu 02 ) 
i? 2 -0 2 -0 2 

+ 2 (i,x)+a 2 -0 2 (i,a;)} . 



(4) 



In the time-independent homogeneous minima of the ac- 
tion or vacua of our model, 4>\ = 4X — 0, 4X = i-R 
(North and South Poles), the Z 2 x Z 2 x Z 2 , a — > 
(— l) <5<,b 06, 6 = 1,2,3 symmetry of the action (fT|) is 
spontaneously broken to: Z2 x Z2, Q — > ( — l) <5a,3 ( 9, 
a, P = 1,2. Finite energy configurations require: 



lim i 

a: — >±oo 







(5) 



Therefore, the configuration space C = {Maps(IR, § 2 )/ 
E < +00} is the union of four disconnected sectors C = 
Cnn U ^ss U £ns U ^sn labeled by the vacua reached by 
each configuration at the two disconnected components 
of the boundary of the real line. 

We now solve the constraint by using spherical coordi- 
nates: 6 G [0,7r], tp G [0,2tt) 

<pi(t,x) = R sin 6(t, x) cos tp(t, x) 
02 (t, x) = R sin 9(t, x) sin tp(t, x) 
03 (t, x) = Rcos6(t, x) . 

In spherical coordinates the mass terms (we shall denote 
in the sequel: a = yl — a 2 ) are 

R 2 



V(6, V) = Y sin 2 0(a 2 + a 2 cos 2 tp) 



(6) 



the action becomes 



S = y dtdx [0„00"0 + sin 2 ed^d^ip] 



R 

— — sin 2 6»((T 2 + d- 2 cos 2 93) 



and the field equations read: 

□ - isin26» (d^ipd^ip - cos 2 99 - cr 2 sin 2 = (7) 

(sin 2 0(9^) - iff 2 sin 2 sin 2^3 = 0.(8) 

Finite energy solutions for which the space-time de- 
pendence is of the form: 

. . . n ( x — vt \ , . ( x — vt 
9(t,x)=0 [-==), p{t,x)=<p 



\/l - it 



for some velocity v, are called solitary waves. Lorentz 
invariance allows us to obtain all the solitary waves in 
our model from solutions of the static field equations 



1 



sin 26 (ip') 



l\2 



1 



cos <p + a sin tp) sin 26 (9) 



-^-(sin 2 ^') = \a 2 sin 2 6»sin2<^ 
dx 2 



(10) 
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where the notation is: 9' = tp' = ^g. The energy of 
the static configurations is: 

E[e,4>] = \ J dx£{e'{x),ip'{x),e{x)^{x)) , 

A/? 2 

£ = — ({9'f + sin 2 9{<p') 2 + sin 2 9(a 2 + a 2 cos 2 p)) 



III. TOPOLOGICAL KINKS 



Equation ^ is satisfied for constant values of <p if and 
only if: tp = 0, S,7r, 4^. Depending on which pair of in- 
constant solution we choose, ([7]) becomes one or another 
sine-Gordon equation: 

□0 + y sin 20 = ; D0+^sin20 = O 

Thus, sine-Gordon models are embedded in our system 
on these two orthogonal meridians. 




Figure 1: a) V(8, tp) deformation of § , b) Embedding of the 
sine-Gordon model at tp = 0, ip = \ as seen in —V(6, tp). 

1. Ki/K* kinks. We denote K 1 /K'l the kink/antikink 
solutions of the sG model embedded inside the S 2 model 
in the tpK x {%) = § or tpx* (%) = two halves of the 
single meridian intersecting the 02 : <^3 plane, 



-IK 



• (x) = 2 arctan e ±ff(x - xo) 



(11) 



see Figure 1. The energy of these kinks, which belong to 
Cns (kinks) or Csn (antikinks), is: E Ki = E K , = 2\R 2 a. 
2. K2/K2 kinks. Taking tpK 2 ( x ) = or tp K *(x) = n > 
we find the sG kinks: 



( x ) = K ,(x) = 2 arctan e ± ( x - xo ' . (12) 



The energy of the K2/K2 kinks, which also belong to 
the Cns, Csn sectors, is greater than the energy of the 
Kx/Kt kinks: E$ 2 = E K » = 2\R 2 . 
3. Degenerate families of Q Q -kinks. When a 2 = 
1, the system enjoys SO(2) internal symmetry and the 
masses of the two pseudo-Nambu-Goldstone bosons are 
equal, there are degenerate families of time-dependent 
Q-kinks of finite energy. If a = 1: ^"(i) = tot + a, 
where ui and a are real constants, solves © for any time- 
independent 9(x). Moreover, by plugging ip® a (t) into ([7]) 
one obtains: 

9 Q * (x) = 2 arctan e ±Vi=^(x-x ) ; V(X ^ 



Therefore, if < lj 2 < 1, the (0 Q ° (x), tp Qa (t)) configura- 
tions form a degenerate circle of periodic in time Q-kink 
solutions of energy: 



E 



2XR 2 



2XR 2 
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Va 



In fact, these Q-kinks can be viewed as the sG kinks ro- 
tating around the main axis of S 2 with constant angular 
velocity to. In another reference frame moving with re- 
spect to the Qo,-kink CM with velocity v, the interplay 
between x and t dependence is more complicated: 



p Qa (x,t) = LU 



t — vx 

VT^v 2 

±0 



a 



9 Wa (x,t) = 2 arctan e 

At the ui = limit we find a circle of static topological 
kinks that form a degenerate family of solitary waves of 
the system. 

Of course, all the multi-soliton, soliton-antisoliton and 
breather solutions of the sG model are also solitons of our 
system in the meridians intersecting either the <p2 ■ </>3 or 
the (pi : cf>3 planes. We shall not discuss these solutions in 
this work and postpone their study to a future research. 



IV. TOPOLOGICAL KINK STABILITY 
A. Small fluctuations on topological kinks 

The analysis of small fluctuations around topological 
kinks requires us to consider both the geodesic deviation 
operator and the Hessian of the potential energy density. 
We will denote 9 = 9 1 6 [0, tt], tp = 9 2 e [0, 2tt), and thus 
the arc-length reads: ds 2 = R 2 d9 1 d9 1 + R 2 sin 2 9 1 d9 2 d9 2 . 
We also denote the Kink trajectories and small deforma- 
tions around them as: 9k(x) = (9 1 K (x) = 9,9 2 K (x) = ip), 
e{x)=e K {x)+r}{x), ti(x) = (r]i(x),rf(x)). 

Let us consider the following contra-variant vector 



fields along the kink trajectory, rj,9' K G T(TS \k)' 
V{x) = r,\x)-^ + r,\x)-^ and 6> K (x) = 0'Jr + tp'^. 

The covariant derivative of rj(x) and the action of the 
curvature tensor on r)(x) are: 



: v = (v H (x) + r) k v j o ,k ) A 



The geodesic deviation operator is: 



R{9' K , v )9' K = 9'W{x)~9' k R\ ]k ^ l . 



D 2 7] 

dx 2 



R(0'k,V)0'k = Ve'Ve-V + R(6'k,v)0'i 



K 



To obtain the differential operator that governs the 
second-order fluctuations around the kink 9k, the re- 
maining ingredient is the Hessian of the potential: 



V^gradV = rj' 



( d 2 V 



r 



dV 



\ d9 l d9i 13 d9 k 



d_ 
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evaluated at 0x(x). In sum, second-order kink fluctua- 
tions are determined by the operator: 

Akv = - (V fl / V e > K ri + R(0'k,v)8'k + V, ; gradU) (14) 



B. The spectrum of small fluctuations around 
Ki/Kt kinks 



Plugging the K\ solutions into (fl"4")) . we obtain the dif- 
ferential operator acting on the second-order fluctuation 
operator around the Kx/K* kinks: 



A Kl r] = A K *r) 



dx 2 



2a 2 



cosh 2 c 



V 



d 



2^,2 



d 2 n 

dx 2 



2atanhax 



drf 
dx 



-2 2 

a rj 



The vector fields v(x) = v 1 (x) 
along the K\ kink orbits satisfy: ; 
, i 



o , 



v 2 (x)(jtfz parallel 



dx 

cotan(2arctane CT 



0, or, 

, v 1 (x) = 1 
, v 2 (x) — coshera; 

Therefore, V\ = , V2{x) = coshera; g|a is a frame 
{i>i,t>2} in r(TS 2 |i<' 1 ) parallel to the K\ kink orbit in 
which (fT5)) reads: 



dx 



d 2 fj 2 

dx 2 



d 2 f) x 



dx 2 



(1 



+ (- 2 

2er 2 
cosh 2 a 



2a 2 



)'/ 



cosh ax 
2 



V 2 



(16) 



where 77 = rj 1 V\ + rj 2 1*2, 77 1 



n 



1 and rj 2 = cosh era; ?7 2 . 



The second-order fluctuation operator (|16[) is a diago- 
nal matrix of transparent Posch- Teller Schrodinger oper- 
ators with very well known spectra. As expected, despite 
the geometric nature of A , we find in the v\ = di- 
rection the Schrodinger operator governing sG kink fluc- 
tuations. Finding another Posch- Teller potential of the 
same type in the 1*2 = -^1 direction comes out as a sur- 
prise because there is no a priori reason for such a be- 
havior in the orthogonal direction. 

In the v\ direction there is a bound state of zero eigen- 
value and a continuous family of positive eigenfunctions: 



ikax / 



sech ax 







f)l(x) = e lkax {t&vhax - ik) 



: W(fc)= C r 2 (fc 2 + 1) 



cosh era; direction the spectrum is similar 



In the V2 

but the bound state corresponds to a positive eigenvalue 



iik( x ) 



7 li-a 2 ( x ) — sech era; 
(tanhcrx — ik) 



ikax 



-(2) _ 



1 - a 2 > 



2;„2 



a 2 k 



Because there are no fluctuations of negative eigenvalue, 
the K\jK\ kinks are stable. 



C. One-loop shift to classical K\/K'{ kink masses 

The reflection coefficient of the scattering waves in the 
potential wells of the Schrodinger operators in (flT))) be- 
ing zero, it is possible to use the Cahill-Comtet-Glaubcr 
formula @ (see also [H for a more detailed derivation) 
to compute the quantum correction to the K\ classical 
kink mass up to one-loop order: 



E Kl (cr) 



E^ i (a) + AE Kl (a) + 0( 1 ^) = 

2XR 2 a [sin v\ H — sin 1/2 — v\ cos v\ 

7r a 



1 



-V2 COS V2\ 



0( 



i? 2 ' 



(17) 



In ()17|) v\ = arccosO = v 2 — arccos cr, are determined 
from the eigenvalues of the bound states and the thresh- 
olds of the continuous spectra. This simple structure of 
the one-loop kink mass shift occurs only for transparent 
potentials. In our model, we find the formula: 



Act 



1 



Ek, (cr) = 2XR a 2 arccos(a) +0(—,) (18) 

7r L a 1 Rr 

For instance, for a = i we obtain a result similar to the 



mass shift of the A0|-kink: 



EkA- 2 ) 



XR 2 - Z ± (I 



2tt V 3 



—) 

6V3/ 



As in the Ac6f-kink case, a zero mode and a bound eigen- 
state of eigenvalue e^' = | contribute. The gaps be- 
tween the bound state eigenvalues and the thresholds 



,(i) 



(0) 



,(2) 



(0) 



1 of the two branches of the 



continuous spectrum are the same in our model. The 
gaps, however, are different from the gaps in the A02 
model between the eigenvalues of the two bound states 
and the threshold of the only branch of the continuous 
spectrum. Both features contribute to the slightly differ- 
ent result. The a = 1 symmetric case is more interesting. 
We find exactly twice the spectrum of the sG kink: two 
zero modes and two gaps with respect to the thresholds 
of the continuous spectrum equal to one. No wonder that 
the one-loop mass shifts of the degenerate kinks is twice 
the one-loop correction of the sG kink: 



Moreover, the quantum fluctuations do not break the 
50(2)-symmetry and our result fits in perfectly well with 
the one-loop shift to the mass of the Af = (2, 2) SUSY 
CP 1 kink computed in [l8| where the authors find twice 
the mass of the Af = 1 SUSY sine-Gordon kink. A dif- 
ferent derivation of formula (|18|) following the procedure 
of see also 0, [2l|, will be published elsewhere. 
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D. The spectrum of small fluctuations around 
K2/K2 kinks 



By inserting the K2 solutions in (|14l) the second-order 
fluctuation operator around the KijK^ kinks is found: 



A K2 i] = A Ki r] 
d 2 V 2 



dx 2 



2 tanh x 





Ki- 


2 


dx 2 


cosh 2 x 


dif 


a 2 V 2 


d 


dx 


09 2 ' 



d 



(19) 



Solving again the parallel transport equations, now along 
the K 2 solutions, it is obtained the parallel frame: 
{ui,u 2 } G T(TS 2 \k 2 ), u\ = gfr, u 2 {x) = coshx^, 
to the K2/K2 orbits. (|T9|) becomes: 



+ 



d 2 f) 2 

dx 2 



d 2 ^ 1 
dx 2 



(1 
2 



cosh x 



cosh x 



I') 1 



(20) 



with 77 = fj Ui + rj 2 U2, T) 1 = fj 1 , rj 2 = cosha^ 2 . 

Again, the second-order fluctuation operator (Til?]) is 
a diagonal matrix of transparent Posch- Teller operators. 
In this case, there is a bound state of zero eigenvalue and 
a continuous family of positive cigenfunctions starting at 
the threshold eW(0) = 1 in the u\ = ~§r direction: 



? 7fcO) 



fjl (x) = sech x 
3 ifca: (tanha; - ik) 



4 X) = 



: {1 \k) = (fc 2 + l) , 



as corresponds to the sG kink. In the U2{x) = cosh x 
direction, the spectrum is similar but the eigenvalue of 
the bound state is negative, whereas the threshold of this 
branch of the continuous spectrum is e^^O) = a 2 : 



Vk( x ) 



Va- 2 -l( x ) = secnx 

= e lkx (tanh x-ik) , e [2 \k) 



e (2) , = a 2 - 1< 
o — 1 



Therefore, K2/K2 kinks are unstable and a Jacobi field 



for k 



10 arises: 



rfj(x) 



"(tanha- — a), = 



E. One-loop shift to classical K2/K2 kink masses 

Once again we use the Cahill-Comtet-Glauber formula 
to compute the quantum correction to the K2 classical 
kink mass up to one-loop order. As before, the angles 
v\ = arccos(O) = -|, V2 — arccos(i<r), are determined 
from the eigenvalues of the bound states and the thresh- 
olds of the continuous spectra. The novelty is that since 
the bound state eigenvalue is negative V2 is purely imag- 
inary. Therefore, we find: 



E K Ja) =2XR 2 - — 

7T 



-a log 



1 IT. o K 

- + \J2 - a 2 - i-a 
a 2 

+ °(^)- ( 21 ) 



The key point is that the one-loop mass shift is a complex 
quantity, the imaginary part telling us about the life- 
time of this resonant state. In the a = 1 symmetric 
case, however, we find the expected purely real answer: 
E K2 (l)=2\(R 2 -l)+0(^). 



F. BPS Q Q -kinks as d = 1 + 1 dyons 

In the fj 2 — 1 case there is symmetry with respect 

to the exp[a( jj ~o jjj] G §0(2) subgroup of the 0(3) 

group. The associated Noether charge distinguishes be- 
tween different Q a -kinks : 

Q = - J dx (4>id t (j)2 ~ <h<h<t>i) = r2 J dx sm2 Odt<P 

Q[Q a ] = R 2 ^ / dx sin 2 9 Qa =2R 2 ^ . 

J u 

For configurations such that 9 is time-independent and 
ip is space- independent, the energy can be written as: 



E 



XR 2 



dx {sin 2 6[<p 



[6' ±<I>sin<9] 2 } 



\R 2 J dx {lu sin 2 9ip =F u8' sin6>} , 



(22) 



(ip = ^-(t)), in such a way that the solutions of the first- 
order equations: 

(p = ui =>• ip Qa (t) = ut + a 



^fu> sm ( 



'(x) = 2arctane =Ft:;(a; - a;o) , 



the Qa _ kinks, saturate the Bogomolny bound and are 
BPS: 



Ebps = — =\{luQ + CoT} 



(23) 



Here, the topological charge T = \ W[6(+oo, 1)] — 
W[0(— oo,f)]| coming from the superpotential W = 
R 2 (l =F cos9) valued at the Q a -kinks gives: T[Q a ] = 
2R 2 ,Va. This explains why "one cannot dent a dyon" 
(even a one-dimensional cousin), see f23j . Conservation 
of the Noether charge forbids the decay of Q a kinks, all 
of them living in the same topological sector, to others 
with less energy. 



G. Bohr-Sommerfeld rule: Q-kink energy and 
charge quantization 

The Bohr-Sommerfeld quantization rule applied to pe- 
riodic in time-classical solutions in our model reads: 



dt I dxw v (t,x) — (t,x) 



R 2 { dt { dx sm 2 9(x,t)^^ = 2nn . 
Jo J dt dt 



6 



In [161 ] it is explained how derivation of this formula with 
respect to the period T = ^ leads to the ODE: A|§ = 

lu-^E), or, 



dn 



EdE 



e VE 2 - 4A 2 i? 4 



E n = X\Jn 2 + 4R 4 



starting from Eq = 2XR 2 and assuming n to be a positive 
integer. The Q-kink energy is thus quantized and the 
frequencies and charges allowed by the Bohr-Sommerfeld 
rule form also a numerable infinite set: 



1 



I7F 



V. THE MASSIVE NON-LINEAR § 2 -SIGMA 
MODEL IN SPHERICAL ELLIPTIC 
COORDINATES 

The secret of this non-linear (l+l)-dimcnsional mas- 
sive § 2 -sigma model is that its analogous mechanical sys- 
tem is Hamilton- Jacobi separable in spherical elliptic co- 
ordinates. This fact will allow us to known explicitly 
not only the kink solutions inherited from the embedded 
sG models, but the complete set of solitary waves of the 
system. 



The spherical elliptic system of orthogonal 
coordinates 



see Figure 2(b). The spherical elliptic coordinates of P 
are half the sum and half the difference of T\ and r 2 : 
2u = n + r 2 , 2v = n - r 2 . u G (R0f,R(-K - f )), 
v G (—R9f, ROf). We remark that this version of elliptic 
coordinates in a sphere is equivalent to using conical co- 
ordinates constrained to S 2 , as defined e. g. in Reference 
12411. We shall use the abbreviated notation: 



. u(t,x) . v(t,x) . 

su = sin — - — , sv = sin — - — , sf = sin j 



R ' R 

sin 2 , sv 2 = sin 2 , sf 2 = sin 2 f 

R R 



and analogously for cu, cv, and cf . To pass from elliptical 
to Cartesian coordinates, or viceversa, one uses: 

<pi{t,x) = — susv , q>3{t,X) = — cucv 
si cf 

<h(t,x) = ±^\A su2 - sf2 )( sf2 - sv2 ) • 

whereas the differential arc-length reads: 

2 2 9 2 

, , su — sv z , 9 su — sv , 9 
dsi 2 = 77 ■ du -\ s • dv 



su 2 - sf 2 



sf z - sv 2 



The spherical elliptic coordinates of the North and South 
Poles, and the foci are respectively: (un,Vn) = (i?#/,0), 
{u s ,v s ) = (R(n - f ),0), (u Fl ,v Fl ) = (R0 f , -R0 f ), 
(u F2 ,v F2 ) ee (ROf, ROf), (u Fi ,v F J ee (R(ir~0 f ),R0f), 
(u^,v P2 ) = (R(n-0f),-RO f ). 



The definition of elliptic coordinates in a sphere is as 
follows: one fixes two arbitrary points (and the pair of 
antipodal points) in § 2 . We choose these points with no 
loss of generality in the form: F\ = (Of,n), F 2 = (0f,0), 
Fi = (tt- f ,0), F 2 = (n- f ,n), f G (0, f ). 

The distance between the two fixed points is d = 2/ = 
2R0f < nR, see Figure 2(a). Given a point Pe§ 2 , let 





Figure 2: a) Foci and antipodal foci of the elliptic system of 
coordinates on § 2 . b) Distances from a point to the foci. 



B. Static field equations and Hamilton-Jacobi 
separability 

We choose a system of spherical elliptic coordinates 
with the foci determined by Of = arccos a, i.e., a 1 = 
cos 2 0f, a 2 = sin 2 0f. We stress that the foci (and their 
antipodal points) are the branching points mentioned in 
the previous Section. In this coordinate system the ac- 
tion for the massive non-linear § 2 -sigma model reads: 



dtdx 



-dnud^u 



-dnvd^v 



V(u,v) 



2(su 2 — sv 2 ) 
The static energy reads: 



su 2 — sf sf — sv 2 

V(u(t,x),v(t,x))} , 

[su 2 (su 2 -sf 2 ) + sv 2 (sf 2 -sv 2 )] 



us consider the distances v\ G [0,7ri?] and r 2 G [0,7ri?] 
from P to F\ and F 2 . 



n = 2R arcsin|/ — (1 — cos#/ cos^ + sinOf sin ^ cos (p) 



E[u,v] = X / dx £ (u'(x), v'(x), u(x), v(x)) , 



sf 



sf — sv' 



-(v'f 



+ V(u, v) . 



r 2 = 2R arcsin y - (1 — cosOf cosO — sinOf sinOcosip) , 



Let us think of E[u,v] as the action for a particle: £ as 
the Lagrangian, x as the time, U(u, v) = — V(u, v) as the 
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mechanical potential energy, and the target manifold 8 2 
as the configuration space. The canonical momenta are: 
P u = Ji7 ' P y = M 7 ' an d the static field equations can be 
thought of as the Newtonian ODE's: 



d 
dx 



su 2 - sf 2 



SV 
Su 
SV 
Sv 



dx V sf ~ — sv 
Because the mechanical energy is 

U{u,v) = -V(u,v) = ^-^ (f(u) + g(v)) = 

su — sv z 

fl 2 [su 2 (su 2 - sf 2 ) + sv 2 (sf 2 - sv 2 )] 

2(su 2 — sv 2 ) 

this mechanical system is a Liouville type I integrable 
system, (see [13 )■ The Hamiltonian and the Hamilton- 
Jacobi equation of spherical Type I Liouville models have 
the form: 



H -- 
dS 



i(su 2 -sf 2 )p 2 -/( U ) 



h v = |(sf 2 -SV 2 )p 2 -g(v) 



dS OS 

ox \ ou ov 



which guarantees HJ separability in this system of co- 
ordinates. The separation ansatz iS(cc, u, v) = —i\X + 
S u {u) + S v (v) reduces the HJ equation to the two sepa- 
rated ODE's, in the usual HJ procedure, leading to the 
complete solution: S = S(x, u, v, i\, i 2 ): 



S = -hx + sgipu) J du^ 



, 2(^+ziSu 2 +/( M )) 
su 2 — sf 2 



w 2( - * d :^; s(v)) m 



the integrations in (|24[) : W(u, v) = S u (u,ii = 0, i 2 — 
0)+S v (v,h=0,i 2 = 0), 

WiPufh) ( u> v ) = F^\u) + G ifi2 \v) 

with (-l)' fll = -sgp„, (-1)' 32 = sgPv ■ sgw, and: 

F^(u) = i? 2 (-l)^cu, G {M (v) = R 2 (-lf 2 cv . 

2. The HJ procedure provides the kink orbits by integrat- 



PV J (sf ^_sv2)|sv| = R3 ^ 2 '- 



tan tan ^7±Z 



tan^| 



tan^| 



tan ~2^~ tan 



(25) 



In Figure 3(a) a Mathematica plot is offered showing sev- 
eral orbits complying with (|2"B"|) for several values of the 
integration constant 72. Note that all the orbits start 
and end at the North Pole and pass through the foci F\ 
such that we have shown a one-parametric family of non- 
topological kink orbits. In fact, there are four families 
of non-topological kinks among the solutions of (j25|) : the 
orbits of a second family also start and end at the North 
Pole but pass through F 2 . The orbits in the second pair 
of NTK families start and end at the South Pole ant pass 
through cither Pi or F 2 . 





K A - f - A F 2 



in terms of the mechanical energy I\ = i\ and a second 
constant of motion: the separation constant I 2 = ■ 



Figure 3: a) Several NTK kink orbits, b) The same NTK 
kink orbits in the elliptic rectangle. 



VI. NON-TOPOLOGICAL KINKS 

We now identify the families of trajectories correspond- 
ing to the values i\ = i 2 — of the two invariants 
in the mechanical system. These orbits are separatri- 
ces between bounded and unbounded motion in phase 
space and become solitary wave solutions in the field- 
theoretical model because the i\ = i 2 = conditions 
force the boundary behavior (See [2|| and [2!| for 
application of this idea to the search for solitary waves in 
other two-scalar field models with analogous mechanical 
systems which are HJ separable in elliptic coordinates.) 
1 . In a first step we find the Hamilton characteristic func- 
tion for zero particle energy (i\ = = i 2 ) by performing 




Figure 4: NTK energy densities for three different values of 
72: 1) 72 = —3, highest peak on the left (blue) 2) 72 = 0, 
symmetrical peaks (green) 3) 72 = 10 highest peak on the 
right (red). 

3. The HJ procedure requires similar integrations in 
S SP" / (si'-tf") _ S SP« / (sf^-sv") = R ( x + 7i) to find 
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the kink profiles (or particle "time " schedules): 



72 are plotted. 



= 2(a;+7i)cf 



, u(x) — f , u(x) + f 

tan 2R tan 2fl 



tan v( - x 2 ) R i tan v 



(26) 4. Reshuffling equations ([25]) and (|26|) . it is possi- 
ble to find the NTK families analytically, (|27|). based 
on (ujv,wjv) = (R9f,0). The other families, based on 
In Figure 4 the NTK energy densities for three values of (us, v$) — (R(tt — 9f), 0) are given by a similar formula. 

I 



tan 



uk (x, 71, 72) 
2R 



±V2y/l + eie 2 t| 



ei + e| + t| +eie2 t 



2 tl 4 - 

2 L 2 



(ei + e 2 2 + t| 4 +e ie 2 t§**)2 - 4(1 + ei) 2 e 2 t 



2 + f 4 \ 



tan 



^A'(a:,7i,72j 
2R 



±J ei+e 2 2 + t| 4 +eie| t| 4 -y(ei + e| + t| 4 +eie 2 ; t| 4 ) 2 - 4(1 + ei) 2 e| t 

V / 2v / T~F^Tti 



2 + f4 
2 



where we have used the new abbreviations: e\ = e 2 ( x +n) c! , e2 = e x+l1 R2 ^ sl2 t^—tan-^. 



(27) 



VII. NON-TOPOLOGICAL KINK 
INSTABILITY: MORSE INDEX THEOREM 

To study the (lack of) stability of NTK kinks, it is 
convenient to use the following notation for the elliptic 



variables: u 
read: 



u, u 



The static field equations 



R.^L = giJ^L (g. k —\ = gV— (28) 

dx dx ^ dx \ jk dx ) ® dui 

Let us consider a one-parametric family of solutions of 
(f28|) : u l K (x;^f). The derivation of 



the subsequent linear system. We skip the (deep) sub- 
tleties of this calculation and merely provide the explicit 
analytical expressions: 



J 



NTK 



(a; 72) 



i? 3 (su 2 -sf 2 )(sf 2 -sv 2 ) 



d d 

sg(Pu) su sg{ Pv ) sv — ) . (29) 

du dv , 



In Figures 5 a)-b), 6 a)-b) two Jacobi fields for two values 
of 72, as well as the corresponding NTK field profiles, are 
plotted for the three <fii, fa, fa original field components. 



D 

dx 



da 



K 



dx 



4) 



dV 



9ik 



du k K 
(?7 



with respect to the parameter 7 implies: 



dx 2 
1 



i:h 



du 3 . 



K 



du k K 



du 



K 



d 2 V 
du 3 du k 



dx dj dx 



R 



jki 



dV 
du 1 



du 3 K 
97 




Figure 5: a) Profiles of the field components for NTK 72 = 
kink, b) Plot of the Jacobi field J NTK (x;0) 



In the last three formulas the metric tensor, the covari- 
ant derivatives, the connection, the curvature tensor, and 
the gradient and Hessian of the potential are valued on 

(m^jU 2 ^), see [HI. Thus, ^j^- is an eigenvector of the 
second order fluctuation operator of zero eigenvalue. The 
derivatives of the NTK solutions (|27|) with respect to the 
parameter 72 are accordingly eigenvectors of the second 
order fluctuation operator of zero eigenvalues orthogonal 
to the NTK orbit, i.e., Jacobi fields that move from one 
NTK kink to another with no cost in energy. 

Better than direct derivation of (f2T)) the Jabobi fields 
can be obtained from (|25p and (f2"6"|) by using implicit 
derivation with respect to the parameter 72 and solving 




Figure 6: a) Profiles of the field components for NTK 72 = 1 
kink, b) Plot of the Jacobi field J NTK (x; 1) 

The zeroes of the Jacobi fields along a given 72- 
NTK orbit (in the four disconnected sectors) are as 
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follows: cither A = (uk{— oo]^/ 2 ) = f)Vk{~ 00572) = 
0), Fi = (ujf (71; 72) = irR - /, 1^(71; 72) = -/), 
•Fa = («Kr(7i;7a) = ttR - f,v K (li;i2) = /), or, 
A = Ju A '(-oo;7 2 ) = ixR - f,v k (-oo;j 2 ) = 0), Fi ee 
(uk(ii,12) = f,v K (^uj 2 ) = -/), F 2 = (uif (71572) = 
/, t'if (71; 72) = /)■ Thus, the conjugate points with re- 
spect to either the North or the South Poles along the 
NTK orbits are listed below: 



Starting Point 


Conjugate Point 


Conjugate Point 


North Pole : A 


Antipodal Focus : F\ 


Antipodal Focus : F2 


South Pole : A 


Focus : Fi 


Focus : Fi 



In this two-dimensional setting, the Morse index the- 
orem states that the number of negative eigenvalues of 
the second order fluctuation operator around a given or- 
bit is equal to the number of conjugate points crossed by 
the orbit (3(j. The reason is that the spectrum of the 
Schrodingcr operator has in this case an eigenfunction 
with as many nodes as the Morse index, the Jacobi field, 
whereas the ground state has no nodes. The Jacobi fields 
of the NTK orbits cross one conjugate point, their Morse 
index is one, and the NTK kinks arc unstable. 



VIII. NON-BPS NON-TOPOLOGICAL KINKS 

The availability of the Hamilton characteristic function 
as a sum of one function of u and another function of v 
allows us to write the energy of static configurations a la 
Bogomolny: 



. . A f, ( su 2 -sv 2 (du su 2 -sf 2 

E\u,v\ = — I dx < tt — ~ ~ — : 

2 J ] su 2 -sf 2 \dx su 2 -su 2 du 



su 2 -sv 2 fdv sf 2 -sv 2 dG^ 



sf — sv 2 V d% su 2 — su 2 dv 

, f , dudF^ , f , dvdG^ 
+ A / dx 1- A / dx - : 

J dx du J dx dv 

Solutions of the first-order equations 
su 2 - sf 2 dF^ 



dv 
dx 



du 



dv sf 2 - sv 2 dG^ 2 



dx 



-Ri-lf 1 

-R(-iy 



su 2 -sf 2 



sf 2 -sv 2 



su(30) 
sv(31) 



are absolute minima of the energy and therefore are 
stable. Note that the energy of the solutions of (|30|) - 

(|3Tj) is positive or zero because sgu' = sg dF ^ l) and 

Sg^ = sg^ 5 ^^. 

Even though the NTK trajectories are solutions of the 
analogous mechanical system provided by the HJ proce- 
dure that is closely related to the ODE system (|3T)|) - (f3"Tj) . 



they do not strictly solve (|30|) - (f3Tj) . Taking the quotient 
of the two equations in (|3"01) - (f3"Tj) we find the equation 



du = gi _ g2 su 2 -sf 2 su 
dv sf 2 - sv 2 sv 



(32) 



which determines the kink orbit flow. Note that this 
equation is identical to the equation in the HJ procedure 
that one must integrate to find l]25p. The subtle point, 
however, is that this flow is undefined, ^, at the four 
foci: Fi, F2, F±, F2, and all the NTK orbits pass through 
one of these dangerous points, see Figures 3(a) and 3(b). 
The non-topological kink orbits solve ([3U|) - PT|) for a given 
sign combination before meeting at a focus and are so- 
lutions of (|3T)|) - (f3"Tj) with another choice of signs after 
leaving these orbit intersections. Thus, non-topological 
kinks arc classified as non-BPS in the terminology of 
"pre-supersymmctric" systems. We remark that in el- 
liptic coordinates the pathology is not in the Hamilton 
characteristic function but in the factors induced by the 
change to elliptic coordinates. The conclusion is that the 
energy of the NTK kinks must be computed piecewise 



along the orbit. 



^(72) 



2A|G^(0)-G(' 32 )(^)| + 



2\\F^\u%) - F^(ug)\ 



££ (72) = 2Ai? 2 |l -a\+ 2Ai? 2 |2cr| = 2Ai? 2 (l + a) (33) 

gives the kink energy as the action of the corresponding 
trajectory. 



A. Singular K\ and K2 kinks: kink mass sum rule 

Analysis of the BPS/non-BPS nature of the topological 
kinks in elliptic coordinates is illuminating. The K\/K* 
kink orbits lie in the v = line, splitting the two-halves 
of the elliptic rectangle: vk x = vk* = 0, sec Figure 3(b). 
The first-order equations (|30|) - (|3ip on the K\jK\ kink 
orbits {fi\ = gives kinks and /3i = 1 anti-kinks) and the 
Ki/K* kink profiles in elliptic coordinates are: 

dx su 



uk 1 (x) = uk'{x) = i?arccos[crtanh((— l) /3l crx)] 
The Ki/Kl kink energy saturates the BPS bound: 

F^\u Kl (+00)) - FW (u Kl (-00)) 



E C K 



A 



2Ai?V 



The K2/K2 kink orbits are the four edges of the elliptic 
rectangle: uk 2 = uk* = R&f, vk 2 = R@fi v k 2 = ~ R&f> 
uk 2 = ukj = — Sf), see again Figure 3(b). The 
K2/K2 kinks are accordingly three-step trajectories in 
the elliptic rectangle. 



I. —00 < x < log tan and u\ 



K 2 



R{-K-6 f ), 



the first-order ODE, and the solutions are: 
j3 2 = 1, v' = R\ sv|, v T K (x) = — v I K .(x) = 2i?arctane 2: 



10 



II. log tan -£ < x < log tan - 2 1 , 



the first-order ODE and the solution are: 

Pi = 0, u' = —Rsu, u 1 ^ Ax) = Uk*(x) — 2i?arctane _a: 



III. log tan — 



< x < +oo, u 



in 



order equation and the solutions arc: 



lift! , the first- 



(3 2 =0,v' = -R\ sv |, vj£*(x) = -v^i{x) = 2i?arctane- a: 

Anti-kinks are obtained by changing the choices of (3\ 
and 02- In any case, the K 2 /K 2 kink energy is not of the 
BPS form: 

E% = A G« (v(-oo)) - G« (v(log tan ^)) 
F(°)(u(logtan^)) - F^Klogcotan^)) 



+A 



+A 



G (0) (u(logcotan-^)) - G (0) (u(+oo)) 



Ai? 2 |l-cf I +XR 2 



2 cf I 



Ai? 2 |l - cf +11 = 2Ai? 2 



It is remarkable that these energies satisfy the following 
"Kink mass sum rule" : 



^(72 



2Ai? 2 (l + cr) = E 



(34) 



In fact, the I72I — + 00 limit of the family of K l2 (NTK) 
kinks is compatible with equation (|25l) only at the edges 
of the elliptic rectangle (forming the K 2 jK^ orbits) and 
the Ki/K* orbit. Therefore, the K x and K 2 form the 
boundary of the moduli space of K l2 in such a way that 
(fM)) shows this combination as one of the NTK kinks. 



A + [<!>(£, x))] is singular on the negative 03-axis but 
a gauge transformation to A~ [&(t, x)]] moves the 
Dirac string -henceforth a gauge artifact- to the 
positive 03-axis. The scalar fields arc constrained 
to live in the <f>' 2 {t,x) + 4>l(t,x) + 4>\{t,x) = R 2 
sphere, a surface where this magnetic flux is con- 
stant. Therefore, Stoke's theorem tells us that Swz = 
R 2 J dx § J2l=i d(j) a (x)A a [$>(x)) is the area bounded by 
a closed curve in § 2 . 

The important point is that the Euler-Lagrange equa- 
tions for the sum of the two actions Swz + S, where S is 
the action of our model, are: 



^ 3 3 



R 



2^ 2^ e *bc<t>c 



6=1 c=l 



dt 



□ c 



+ f|0c 



(35) 



At the long wavelength limit, the ODE system (|35|) be- 
come the Landau-Lifshitz system of equations of ferro- 
magnetism. The connection between the semi-classical 
(high-spin) limit of the Heisenberg model and the quan- 
tum non- linear § 2 -sigma model is well established |32j. 



A. Spin waves 

Plugging the constraint into (f3"5)) . we find the system 
of two ODE's: 



Sg03 



IX. SOLITARY SPIN WAVES 

Field configurations that satisfy the Euler-Lagrange 
equations: 

3 



dAa 

dt 



3 3 

=EE^(t.«)]-j(^) 

6=1 c=l 

are extremals of the "Wess-Zumino" action: 
Swz[$] = R 2 I dtdx MHKx)]^{t, 



a=l 



In particular a "magnetic monopole " field B a [$(t,x)] = 
^"^j' 3 ^ in the M 3 internal space where the § 2 -sphere is 
embedded is obtained by the choice of singular "vector 
potentials" : 

Afm, X )} - 



y/<fi + <% + <%(<h ± y/<% + 01 + 0|) 



V# + $ + ± V01 + 02 + 0!) 



A*[$(t,x)] = 



+0o 



E 7 07 d < 



R 2 -J2 



+ D0/5 + "^0/3 



7 V7V7 



0/3 



i? 2 -E 



7 



E 7 7 9 M 7 + Ea 0.5^05 



i? 2 -E 



^ (9 /i 7 a M 7 + 7 n0 7 ) 



7 <^7^7 



. 



(36) 



ce,/3, 7 = 1,2, m 2 = 1, = a . The ground states are 
the homogeneous solutions of this system: 0° = 0" = 0, 
0§ = ±i?. In order to visualize these configurations in, 
c. g., Figure 7 we draw the spin chain in such a way that 
the 02 : 03 plane is perpendicular to the x spatial line 
whereas 0i is aligned with the x-axis. We stress that this 
choice of basis is arbitrary but it is easy to figure out the 
formulas and the graphics in another rotated basis for 
the magnetization vector: <p(x) = 4>i(x)e\ + <f> 2 {x)e 2 + 
03(a^)e3 = (f>i(x)e'i + 2 (:r)e / 2 + 3 (a;)e'3. The main fea- 
tures of our preferred basis ei, e 2 , e*3 are: 1) The e\ vector 
points in the direction of weaker —V{6, (p) potential, see 
Figure 1(b). 2) ei,e 2 ,e 3 is the basis used in the contin- 
uous XY (in fact YZ) model of easy-axis ferromagncts 
near the Curie point, see [H, [3f| and References quoted 
therein. 
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Figure 7: a) Ground state <f>% = R . All the spins are aligned 
pointing to the North Pole b) Ground state 4>% = —R. All the 
spins are aligned pointing to the South Pole. 




Figure 8: Graphic arrow representation of the K\ kinks: a) 
K\ spin chain, b) Perspective from one component of the 



The spin fluctuations 4>i(t, x) = 8<pi(t, x), 4>2{t,x) = 
5({>2(t,x) around the ground state cj)z{t,x) = R satisfy 
the linearized equations: 







dSch , d 2 Scj) 1 d 2 6(j) 1 
dt 







dt 2 

dt + dt 2 
Therefore, the spin waves: 



3a; 2 
d 2 S> 



dx 2 



b a (t,x) = — L= V — 7^== (a a (k)t 
V XL ^ yfu{k) 

+a* a {k)e- lut+tkx ) 



iut—ikx 



(37) 



satisfying periodic boundary conditions S<f) a (t,x) = 
Scf> a (t,x + XL) are solutions of |37|) for the frequencies 
complying with the homogeneous system of algebraic 
equations: 



-LU 2 + k 2 

—iuj 



1 ioj 

-uj 2 + k 2 + a 2 




At the long wavelength limit uj 2 << u>, (|38|) is tanta- 
mount to the non-relativistic dispersion law 

w 2 (k 2 ) = (k 2 + l)(k 2 + a 2 ) 

characteristic of ferromagnetic materials, although the 
quadratic terms in the free energy prevent the standard 
u)(k) = k 2 form. 



means of a 7r-rotation around the 0i-axis. 



and are also spin solitary waves in this system, (Fig. 9). 




Figure 9: Graphic arrow representation of the K2 kink a) 
K2 spin chain, b) Perspective from one component of the 
boundary of § 2 x R showing a forward spin flip. 



Because the system of ODE's giving static solutions 
of the ([33)) PDE system is the same as the static field 
equations of the non-linear S 2 -sigma model, the NTK 
kinks are also solitary spin waves, see Figure 10. 




B. Bloch and Ising walls 

One may check that the K\jK\ kinks (fTTj) solve the 
static Landau-Lifshitz equations (|36|) on the <f>\ = orbit: 



<i 2 <. 



dx 2 



i? 2 -0 2 



d<t>2 \2 



R 2 



dx 



dx 2 



The K\ /K"l kinks of the non- linear sigma model are con- 
sequently solitary spin waves of this non-relativistic sys- 
tem, see Figure 8. 

Simili modo, the K^jK^ kinks (JT^J) solve (f3"tj|) along 
the 60 — kink orbit: 



dx 2 R 2 



dx I 



R 2 



dx 







1 dx 2 



Figure 10: Graphic arrow representation of K 12 kinks: a) 
K 12 spin chain, b) Perspective from the boundary of S 2 x R 



showing the 2ty rotation around the 
come back to the initial ground state. 



>i-axis of the spin to 



In sum, understood as solitary spin waves K\jK\ kinks 
arc Bloch walls whereas K^jK^ kinks are Ising walls de- 
scribing interfaces between ferromagnetic domains, see 
[33l ] . [3a |. In this model we have thus found a mod- 
uli space of solitary waves with an structure very sim- 
ilar to the structure of the space of solitary waves of the 
XY model described in References [33| and [HI- There 
are Bloch and Ising walls and a one-parametric family 
of NTK kinks that are non-linear superpositions of one 
Bloch and one Ising wall with arbitrary separation be- 
tween their centers. The novelties here are: a) there is 
no need in the free energy of fourth-order terms in the 
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magnetization in the non-linear sigma model for finding 
these mixtures of Bloch and Ising walls, b) The analyti- 
cal expressions (|27|) differ from their analogues in the XY 
model. 

From the stability analysis performed in previous Sec- 
tions, it is clear that only the Bloch walls are stable and 
saturate the Bogomolny bound. Things are different at 
the (7=1 limit where all the kinks are topological, Bloch 
walls, and saturate the Bogomolny bound. In this lat- 
ter case the structure of the kink space is akin to the 
kink space structure of the BNRT model [IB], see [37l ]. 
[38| , [HI ■ There is a one-parametric family of degenerate 
Bloch walls saturating the Bogomolny bound. 



X. FURTHER COMMENTS: 
SUPERSYMMETRY AND STABILITY 



Finally, we briefly explore the possibility of embed- 
ding our bosonic model with its moduli space of kinks 
in a broader supcrsymmetric framework, ft turns out 
that the simpler Af = 1, d = 1 + 1 SUSY version of 
the massive non-linear § 2 -sigma model only exists if the 
masses of the pseudo Nambu-Goldstone bosons are equal 
((j=l). It also seems difficult to build more exotic pos- 
sibilities coming from dimensional reduction of models of 
Kahler or hyper-Kahler nature because the potential en- 
ergy density is not compatible with complex structures 
when a =/= 1. 



and we rewrite the second order fluctuation operator 
around the K\ kink (with \k 1 ( x ) = Re~ ax )) in the form: 

A^ = -^ + 2.(l-tanha,)^ 

- (1 - 2a 2 + 2a 2 tanhax) r/) A - (^L + 

drj 2 _ 2 ^ oa 1 N 2 \ d 



-2a {1 — tanhera;) — h a" (1 — 2tanh<72i) rj , 

dx J ox 



In a parallel frame fi = fj}{x)-^r + /i 2 (x)^p- € 
r(T§ 2 \ Kl ) , ^ + Y) k { X K)X%V k = 0, along the K 1 kink: 

^+a(l-tanh)fj, 1 (x) =0 =S> ^(x) = 1 + e~ 2ax 
dx 

$L+a(l- tanh)/i 2 (s) = => ^(x) = 1 + e~ 2ax . 
dx 



we recover the Posch- Teller operators: 



d 2 V 2 
lb? 



dx 2 
+ (- 2 



2cr 2 , ,\ , „„. d 

1 rS—W) l + e- 2ra hrr 

cosh era; / ox 

tf) (1 + ^)^1 (41) 



cosh 2 



Note that now the if i orbits are the positive and negative 
ordinate half-axes, the stereographic projections of the 
ip = § and <p = 4p half-meridians, such that fluctuations 
orthogonal to the orbit run in the direction of the abscissa 
axis. 



A. Isothermal coordinates 



B. The N — 1 massive SUSY sigma model 



It is convenient to introduce isothermal coordinates in 
the chart S 2 — {(0, 0, —R)}, which are obtained via stere- 
ographic projection from the South Pole: 



R<t>i 



x 



l + f R + &g{fa)y/ IP 

4>2 _ R<p2 

1 l+f " R + SgfaW 



(39) 



The metric and the action in this coordinate system read: 
4R 4 



ds 2 = 



{R 2 +X 1 X 1 +x 2 x 2 ) 2 



(dx'dx 1 +d X 2 d X 2 ) 



2R 4 



S[x ,x] ~J dx (R 2 + x 1 x 1 + x 2 x 2 ) 2 
+^xW-(xV + ^Vx 2 )] , 



whereas the K\ kinks are given by: 

Xk, (s) = , X 2 Kl [?) = ±R cxp[±a(a; -*„)], (40) 



In Reference [40j we analyzed the relationship of the 
complete solution of the Hamilton- Jacobi equation for 
zero energy and the superpotential of a supersymetric 
associated classical mechanical system. Thus, we are 
tempted to use the Hamilton characteristic function 

W R 2 + x 1 x 1 + x 2 x 2 _ 

y/ (a + (/3 2 )R 2 + a_(/3 2 )( X 1 X 1 + X 2 X 2 )f - 4^W, 

o-±(/3a) = 1 ± (-l) A cr, to build the TV = 1 SUSY exten- 
sion of our massive non-linear S 2 -sigma model. On one 
hand we have that: 

1 j-dW^'M d W U>i,h) _ 2B 2 (x 1 x l +a 2 X 2 X 2 ) 
2 9 d X l ' dx> ~ {R 2 + XV + X 2 X 2 ) 2 ' 

V/3i,/?2- On the other hand (|32"|) is free of branch points 
only for a = 1. Supersymmetry does not allow superpo- 
tcntials with branch points and it seems that Hamilton- 
Jacobi characteristic functions are compatible with a 
weaker form called pseudo-supersymmetry in [4lj . We 
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close our eyes to this fact for a moment and proceed to 
formally build the AT = 1 SUSY extension of our model 
using ([t2j) . 

There are also two Majorana spinor fields: 

* V) -te) ,( * ) *"*'°" u ' 

We choose the Majorana representation 7 = a 2 ,^ 1 — 
ia\j 5 = a 3 of the Clifford algebra {7^,7^} = 2g» v and 
define the Majorana adjoints as: ip l = ('0 I )*7°. The 
action of the supersymmctric model is: 



S 



dx 2 
2 

! D 7i/7 7i ,fc adWdW Tl DdW ., 
6 OX ax- 7 <J\ ox 3 

where = 357^ + f^r • Thc spinor supercharge 

d*fl« (fW^X J ' + *7W*S (43) 

acts on the configuration space and leaves the action in- 
variant. Time- independent finite energy configurations 
complying with 



da; 



^x 7 



(44) 



annihilates the supercharge combination Qi + Q2 and 
these solutions might be interpreted as i BPS states in 
this supersymmctric framework. In particular, the SUSY 
K\ kinks 



-1 



satisfy 1(33]) (with appropriate choices of /3i, /i^)- Note 
that VkiC 23 ) i s the SUSY partner of Xk 1 ( x ) under the 
action of the broken SUSY supercharge Qi — Q2. We 
also remark that 



d x 2 Kl 



±aRe ±ax = ±aR{l + e ±2ax ) ■ 



1 



dx ' cosh era; ' 

i.e., the fermionic partner in the SUSY kink is the zero 
mode of the second order fluctuation operator back from 
the parallel frame to the K\ orbit. 

C. Fermionic fluctuations 

The Dirac equation ruling the small fermionic fluctua- 
tions on the K\ kink reads: 

D8i/>%x) = i(7°a -7 1 5i)^(i,a;) 

DdW 

+9 %J (x Kl )g^(XK^ k (t,x) (.45) 



Acting on (|45|) with thc adjoint Dirac operator, the 
search for solutions of D^DSip 1 ^, x) = of the stationary 
form 5ip l (t, x) = e lwt Sg l (x, to) requires us to deal with the 
following ODE system: 



it 

dx 



2 5g\x) 



- DdW 

dx 3 dx k 



+ R jkl-^X j Kl -^x'kM(x) 



-»7 



dW 



dx' 



M 



D 2 dW 



dx 3 dx k dx l dx m 



Sg m (x) = u 2 Sg i (x) 



valued at x = Xk x ■ 

On eigenspinors of — ij 1 = a , 5g\ (x) = ±Sg l 2 (x) = 
5g l ± (x), the above spectral ODE system reduce to the 
(symbolically written) pair of equations: 



At 8q± 



d 2 

■-To + W " ® W " 

dx z 



R±W' ® W" 



(46) 

is exactly equal to the second order differential op- 
erator ruling the bosonic fluctuations. Therefore, in the 
parallel frame to the K\ orbit we write A^- in matrix 
form: 



.JL 

dx 2 



1 - 




cosh era: 





v 2 



dx- 



cosh-cra: 



In the same frame A K is the intertwined partner, see 
El: 



A 



.JL + 1 

dx 2 + 1 







d A 

dx- 



If a 7^ 1, there is a bound state in A^ of energy 
1 — g 2 unpaired with an eigenstate of the same energy 
in A^ i , a fact incompatible with supersymmetry as we 
expected from the use of the complete solution of the 
Hamilton- Jacobi equation as superpotential, closing our 
eyes to the fact that, related to the instability of NTK 
and K2 kinks, the Hamilton characteristic function has 
branching points at the foci defining the elliptic coor- 
dinate system. A similar problem arouse in [42j and 
[43| where meromorphic Hamilton characteristic func- 
tions have been found. It is an open problem to explore 
whether or not these milder singularities allow the use of 
these Hamilton characteristic functions as superpoten- 
tials to extend the bosonic models dealt with in [43] , [IH 
to the supersymmetric framework. 

If the masses are equal (c = 1), however, the Hamilton 
characteristic function is free of branching points and thc 
unpaired states are zero modes. Thc M = 1 SUSY model 
is correct and we can apply the SUSY version of the 
Cahill-Comtet-Glauber formula proposed in [44[ to find 
the same one-loop correction to the SUSY S 2 kink as 
given in [l8| : 

A< s > = 1) = — B sin ^ - v t cotf) = -- ■ 



1=1 
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Here i/y = = arccos(O) = -| arc the angles obtained 

from the bound states of At- . There are no bound states 

Hi 

in the spectrum of A^- . 
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